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Origin of the Large Perturbative Corrections to Higgs Production at Hadron Colliders 
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The very large _R'-factor for Higgs-boson production at hadron colliders is shown to result from 
enhanced perturbative corrections of the form {Cat:ols)^ , which arise in the analytic continuation 
of the gluon form factor to time-like momentum transfer. These terms are resummed to all orders 
in perturbation theory using the renormalization group. After the resummation, the /('-factor for 
the production of a light Higgs boson at the LHC is reduced to a value close to 1.3. 



I. INTRODUCTION 

The discovery of the Higgs boson is the most important 
goal of modern particle physics. The inclusive production 
cross sections for pp ^ H + X and pp ^ H + X have been 
calculated a long time ago at next-to-leading order (NLO) 
in perturbation theory P, Q , and since a few years results 
at next-to-next-to-leading order (NNLO) are available [1, 
0, Q . Inclusive Higgs-boson production is thus one of the 
best studied processes from a theoretical perspective. 

In view of this fact, it is uncomfortable that the behav- 
ior of QCD perturbation theory appears to be rather poor 
in this case. The X-factor for Higgs-boson production, 
defined as the prediction for the cross section normalized 
to the Born approximation, takes surprisingly large values. 
For the production of a hght Higgs boson {niH < 150 GeV) 
at the LHC, one typically finds K Ki 1.7-1.9 at NLO and 
K « 2.0-2.2 at NNLO. Also, the residual dependence on 
the renormalization and factorization scales remains signif- 
icant even at NNLO. The standard argument that the large 
JsT-factor results from the accessibility of new production 
channels beyond the leading order, such as qg — > Hq and 
qq — > H, does not apply in this case, as these contributions 
to the cross section are known to be below 10%. Also, the 
/^-factor is not much reduced by soft-gluon resummation 
near the partonic threshold . 

In this Letter we show that the bulk of the large per- 
turbative corrections to Higgs-boson production via gluon- 
gluon fusion originate from terms of the form (C^Tras)" 
arising from the analytic continuation of the gluon form 
factor to time-like momentum transfer, and that these 
terms exponentiate to leading order. 



II. TIME-LIKE GLUON FORM FACTOR 

The Higgs-boson production cross section at hadron col- 
liders such as the Tevatron or LHC is dominated by the 
gluon-giuon fusion process gg ^ H via a top-quark loop. 
For a not too heavy Higgs boson, this process is well ap- 
proximated by the effective local interaction Q 



where v ~ 246 GeV is the Higgs vacuum expectation value, 
and the short-distance coefficient Cj = as/{l2n) + 0{al) 
is known to NNLO Q and has a well behaved perturba- 
tive expansion for fj. mn. The production cross section 
is related to the discontinuity of the product of two such 
effective vertices. It can be written as the convolution of a 
hard-scattering kernel with parton distribution functions. 

The large corrections we identify are due to virtual cor- 
rections to the effective ggH interaction ^ and arise from 
quantum corrections characterized by the scale fj, ~ mn- 
These effects are described by a universal factor and af- 
fect differential distributions in same way as the total 
cross section. They can be factorized into a hard function 
H{rrijj, /i^), which is the square of the on-shell gluon form 
factor evaluated at time- like momentum transfer q^ = mjj, 
and with infrared divergences subtracted using the MS 
scheme d, [13, [HI- On a technical level, the hard func- 
tion appears as a Wilson coefficient in the matching of the 
two-gluon operator in ([T]) onto an operator in soft-coUinear 
effective theory (SCET) [13, [H, in which all hard modes 
have been integrated out. This matching takes the form 



(2) 



where = —q^ is (minus) the square of the total momen- 
tum carried by the operator. The fields .A^^ and A'^'^ are 
effective, gauge-invariant gluon fields in SCET [3|. They 
describe gluons propagating along the two light-like direc- 
tions n, n defined by the colliding hadrons. 

The two-loop expression for the Wilson coefficient Cs 
can be extracted from the results of [l^j. We write 



An 



where L = ln((5^//i^). The one-loop coefficient reads 



ci(L) = Ca{-L 
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(3) 



(4) 



and the result for the two-loop coefficient can be found in 
[lOl [l6| . The hard function is given by the absolute square 
of the Wilson coefficient at time-like momentum transfer, 
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The Wilson coefScient obeys an evolution equation, which 
reflects the renormalization properties of the effective two- 
gluon operator in SCET. It reads [§| 



dlnfi 



- cusp 



(a,) ln^+7^(a,) 
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where F^^^p is the cusp anomalous dimension of Wilson 
lines with light-like segments in the adjoint representation 
of SU{Nc). It controls the leading Sudakov double loga- 
rithms contained in Cs and is known to three-loop order 
The single- logarithmic evolution is controlled by the 
anomalous dimension 7'^, which can be extracted from the 
infrared divergences of the on-shell form factor 9]. Us- 
ing results from 18] it can be derived to three-loop order 
[la |. The evolution equation © links the coefficients of 
the logarithmic terms in ([3]) to coefficients in the perturba- 
tive expansions of the anomalous dimensions and the QCD 
/3-function. At one-loop order we have 



(7) 



where Tq = 4Ca and 7o = 0. 

The Wilson coefficient at space-like momentum transfer 
has a well behaved expansion in powers of the coupling 
constant, if the renormalization scale is taken to be of order 
the natural scale, /it^ ~ Q^. For instance, with Nc — 3 
colors and Uf = 5 light quark flavors, we find 



CsiQ^, Q') = 1 + 0.393 a,(Q2) _ 0.152 al{Q^) 



(8) 



The nature of the expansion changes drastically when the 
same coefficient is evaluated at time-like momentum trans- 
fer = —q^ — ie. We then obtain 

Csi~q^,q^) = 1 + 2.75 + (4.84 + 2.07z) a^q^) 

+ .... (9) 

The expansion coefficients are more than an order of mag- 
nitude larger than in the space-like region. The origin of 
this effect is that the Sudakov (double) logarithms con- 
tained in the coefficients c„(L) in ^ give rise to tt^ terms 
when we analytically continue L ln.{q^ / fi^) — in. For the 
hard function entering the Higgs-boson production cross 
section, this implies 

H(m]j,m]i) = 1 + b.h()as{m]j) + 17.24a^(m|^) -|- . . . 

= 1 + 0.623 + 0.221 + ..., (10) 

where the numerical estimates in the last line refer to the 
NLO and NNLO corrections for a Higgs-boson mass of 
120 GeV, and we use as{m?^) = 0.118 as our normalization 
of the running coupling constant. These hard matching 
corrections account for the bulk of the X-factors found at 
NLO and NNLO. 

The large expansion coefficients in the perturbative se- 
ries for the Wilson coefficient in the time-like region can be 
avoided if we evaluate this coefficient at a time-like renor- 
malization point, in which case (here and below, negative 



arguments of the running coupling are always understood 
with a —ie prescription) 
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(11) 



with L = ln((7^//i^) and the same expansion coefficients as 
in We then obtain 

Cs(-g^ -q^) = 1 + 0.393a,(-q2) _ 0.152 al{-q^) + ... 

(12) 

instead of The perturbative series analogous to that 
in pI7|) reads 

-m|^)|2 = 1 + 0.0845 - 0.0015 + . . . , (13) 

which indeed exhibits a vastly better behavior. 

In the expressions above, the running coupling is evalu- 
ated at time-like momentum transfer — /i^ — ie. The func- 
tion ^^(/i^) in perturbation theory is analytic in the com- 
plex /i^ plane with a (physical) cut on the negative real 
axis and a (unphysical) Landau pole at /i^ = A|^. Since 
we are interested in very large values, the Landau pole 
is not of concern to our discussion. The definition 



I3{as) = 2 



das(pL^) 



n=0 



din /i^ 

of the QCD /3-function implies that 

da _ in 



00 



/3(a) 



(15) 



and this relation allows us to define the running coupling 
at time-like argument in terms of that at space-like mo- 
mentum transfer. At NLO we obtain 



as(-M^) 



l-m(/i^) + ^ 



/3o 4^ 



In [I - ia{^^'')\+0{al) , 
(16) 

where a(/i^) = /3Qas(/i^)/4. In standard applications of 
the renormalization group one would count this quantity 
as an 0{\) parameter. Since numerically a{rnj^) « 0.2, it 
is however also reasonable to count a — 0{as). 



III. RESUMMATION 

What is needed for the calculation of the Higgs-boson 
production cross section is the Wilson coefficient at posi- 
tive, not negative fj,^, see ([5]). We will use the solution to 
the renormalization-group equation ([6]) to relate this coef- 
ficient to the one in pT|) . In that way the large corrections 
arising in the time-like region are resummed to all orders 
in perturbation theory. We write the solution in the form 

Himl,^l^)^U{ml,^l^)\Csi-mj„-^,^)\\ (17) 

where [T9| 



lnU{m'jj,ij'^) = 2Re 



2S{-fi ,Ai ) - a^s{-fj. ,/u ) 
2 



(18) 
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and similarly for the function a^s. The perturba- 
tive expansions of these functions obtained at NNLO in 
renormalization-group improved perturbation theory can 
be found in [2^. They can be simplified using relation 
p6)) . To leading order we find 



In U{mH, fj?) = — < — - — [2a arctan(a) - ln(l + a^)] 



^ 



/3o 



^ 



ln(l + a2) 



(20) 



+ ^ karctan^fa) - In^fl + a^)] + 0{a,) \ , 
4po J 

where a = a(m|^). Note that the result is ^-independent at 
this order. The relevant anomalous-dimension coefficients 
are Fq = 4Ca, 7o = 0, and 



67 TT' 
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(21) 



where Ca — Nc, Tp ~ 1/2, and = 5 is the number 
of light quark flavors. The coefficients of the /3-function 
follow from (HI]) . 

The expression for the evolution function simplifies con- 
siderably if we treat a{m?fj) w 0.2 as a parameter of order 
Us- Inserting the values of the one- loop anomalous dimen- 
sions from above, we then find 



In U (m|j, /i^ 



1 



pA 
^ 



4tt 



0{al 



(22) 

This result makes explicit that the "7r^-enhanced" correc- 
tions are terms of the form (C^Tras)" in perturbation the- 
ory and exponentiate at leading order. The simplest way 
to implement our resummation in existing codes for Higgs- 
boson production would be to multiply the fixed-order re- 
sult with exp[Cyi7ras(m|^)/2] and subtract the expanded 
form of this factor from the perturbative series. This treat- 
ment is sufficient for practical purposes. 

Numerically, setting fj, = ran = 120 GeV we obtain 
InC/ = {0.563,0.565,0.565} at LO, NLO, and NNLO from 
the exact expression for the evolution function derived from 
psp . indicating that the leading-order terms give by far 
the dominant effect after renormalization-group improve- 
ment. The analytical expressions (PO)) and provide 
accurate approximations to the exact results. The first 
equation gives lnJ7 = 0.562, while the second one yields 
lnJ7 = 0.567. The close agreement of these two numbers 
shows that the running of coupling constant between 
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FIG. 1: LO (light), NLO (medium), and NNLO (dark) pre- 
dictions for the Higgs-production cross section at the LHC in 
fixed-order perturbation theory (left) and after resummation of 
the TT'^-enhanced terms (right). 

and — /i^ is a minor effect compared with the evolution 
driven by the anomalous dimension of the effective two- 
gluon operator in ([2]). 

We are now in a position to discuss our improved results 
for the hard function in the formula for the Higgs-boson 
production cross section. Setting fi = ruH = 120 GeV, we 
obtain 

H{m]j,'m]j) = {1.756 (LO), 1-907 (nlo), 1-906 (nnlo)} - 

(23) 

This should be compared with the poorly converging series 
H — {1, 1.623, 1.844} obtained using fixed-order perturba- 
tion theory. Figure[l]illustrates the impact of the resumma- 
tion of the TT^-enhanced terms on the cross-section predic- 
tions for Higgs-boson production at the LHC. The bands in 
each plot show results obtained at LO, NLO, and NNLO 
using MRST2004 parton distributions [2l|. Their width 
refiects the scale variation obtained by varying the factor- 
ization and renormalization scales between m///2 and 2m h 
(setting /ir = ^/). The convergence of the expansion and 
the residual scale dependence at NLO and NNLO are much 
improved by the resummation. The new LO and NLO 
bands almost coincide with the NLO and NNLO bands in 
fixed-order perturbation theory, and the new NNLO band 
is now fully contained inside the NLO band. 



IV. DRELL-YAN PRODUCTION 

The cross section for the Drell-Yan process receives the 
same type of 7r^-enhanced corrections as the Higgs-boson 
production cross section, however in this case no anoma- 
lously large ii'-factors arise at NLO and NNLO. Let us 
briefly discuss why this is the case. 

The vector-current matching coefficient Cy appearing in 
the Drell-Yan case is defined in analogy with Cs in ([2]), but 
with the two-gluon operator replaced by the electromag- 
netic current q^^q [ol. llCt [T]|. It obeys an evolution equa- 
tion of the same structure as ([S]), in which the cusp anoma- 
lous dimension in the adjoint representation is replaced by 



4 



that in the fundamental representation, and in which an 
anomalous dimension 7^ replaces 7'^. The cusp anomalous 
dimensions are simply related by T^gp/F^gp = Cf/Ca 
[13]. The one-loop coefHcient of 7^ is 7^ = -6Cj.. The 
resummation of 7r^-enhanced terms is accomplished as in 
(fT7| . where in the explicit forms of the evolution function 
the appropriate expansion coefficients must be used. 

The one-loop matching contribution to the Wilson coef- 
ficient Cv reads 

cr(i) = -^f i^-f i + C^(y-8), (24) 

which has the same structure as However, the lead- 
ing TT^-enhanced terms in Drell-Yan production are smaller 
than those in Higgs-boson production by a factor 4/9, and 
the constant —8 in the matching condition ((24|) . which is 
absent in the Higgs case, has the effect of partially com- 
pensating the large, positive tt^ terms resulting from the 
analytic continuation L ln(g^/ n^) — itt. This can also be 
seen by looking at some numerical values. For the squared 
coefficients at time-like momentum transfer, we obtain for 
q ~ 120 GeV the perturbative expansions 

|Cy(-g^g2)P = 1 + 0.0845 + 0.0292 +... , 

(25) 

|Cv(-g^ -(7^)1^ = 1 - 0.1451 - 0.0012 +... , 

which should be compared with (fTO| and (|13p . While the 
convergence is better for fi^ — —q^, the large one-loop cor- 
rection seen in the Higgs case is absent for the reasons 
mentioned above. 

It was shown in {22] that the bulk of the perturbative 
corrections to the Drell-Yan cross section are related to tt^- 
enhanced terms in the ratio of the time-like to space-like 
Sudakov form factors. In this paper a resummation formula 



was derived for the ratio of these form factors, which shares 
some similarities with our results (|20|) and ([22|) . To the best 
of our knowledge, a corresponding analysis has not been 
performed for the Higgs-boson production cross section. 
We expect that the approach based on effective field theory 
presented here can be adapted straightforwardly to other 
processes. 

V. CONCLUSIONS 

We have shown that the large A'-factor for Higgs-boson 
production at hadron colliders results from a simple kine- 
matical effect: the analytic continuation of the gluon form 
factor to time-like momentum transfer. This leads to large 
perturbative corrections of order (CATraj)", which can be 
resummed to all orders by solving a renormalization-group 
equation. Our approach employs methods from effective 
field theory and properties of the QCD running coupling 
in the complex momentum plane. 

After the resummation of the 7r^-enhanced terms, the K- 
factor for Higgs-boson production at the LHC is reduced 
to a rather modest value of about 1.3 at both NLO and 
NNLO. A detailed analysis of the phcnomcnological conse- 
quences of our observation, combined with state-of-the-art 
results for NNLO corrections and soft-gluon resummation, 
will be presented in a forthcoming paper fl^. Extensions 
of our approach to other hard-scattering processes with 
time-like momentum transfer, such as event shapes or tt 
production at hadron colliders, will be discussed elsewhere. 

Acknowledgments: We are grateful to Babis Anastasiou, 
Martin Beneke, Frank Petriello and Giulia Zanderighi for 
useful discussions. The research of T.B. was supported by 
the U.S. Department of Energy under Grant DE-AC02- 
76CH03000. Fermilab is operated by the Fermi Research 
Alliance under contract with the Department of Energy. 



S. Dawson, Nucl. Phys. B 359, 283 (1991). [13 
A. Djouadi, M. Spira and P. M. Zerwas, Phys. Lett. B 264, 
440 (1991). 

R. V. Harlander and W. B. Kilgore, Phys. Rev. Lett. 88, [14 
201801 (2002) arXiv:hep-ph/0201206 . 
C. Anastasiou and K. Melnikov, Nucl. Phys. B 646, 220 [15 
(2002) arXiv:hep-ph/0207004 . 

V. Ravindran, J. Smith and W. L. van Neerven, Nucl. Phys. [16 
B 665, 325 (2003) arXiv:hep-ph/0302135 . 

S. Catani, D. de Florian, M. Grazzini and P. Nason, JHEP [17 
0307, 028 (2003) ar Xiv:hep-ph/03062 11 . 

T. Inami, T. Kubota and Y. Okada, Z. Phys. C 18, 69 [18 
(1983). 

K. G. Chetyrkin, B. A. Kniehl and. M. Steinhauser, Phys. [19 
Rev. Lett. 79, 353 (1997) aFXiv:hep-ph/9705240 . 
T. Bec her and M . Neubert, Phys. Rev. Lett. 97, 082001 [20 
(2006) [arXiv:hep-ph /0605050 . ' 

A. Idil bi, X. d. Ji and F. Yuan, Nucl. Phys. B 753, 42 [21 
(2006) [arXiv:hep-ph706 05068 . 

T. Becher, M. Neubert and G. Xu, ,]HEP 0807, 030 (2008) 
arXiv:0710.0680 [hep-ph]]. [22; 
C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. Rev. D 
65, 054022 (2002) arXiv:hep-ph/0109045. . 



C. W. Bauer, S. Fleming, D. Pirjol, I. Z. Rothstein 
and I. W. Stewart, P hys. Rev. D 66, 014017 (2002) 
arXiv:hep-ph/0202088 . 

R. J. Hill and M. Neubert, Nucl. Phys. B 657, 229 (2003) 
arXiv:hep-ph/021 1 018 . 

R. V. H arlander, Phys. Lett. B 492, 74 (2000) 
a rXiv:he p-ph /0007289 . 

V. Ahrens, T. Becher, M. Neubert and L. L. Yang, in prepa- 
ration. 



691, 129 (2004) 


arXiv:hep-ph/0404111_. 


S. Moch, J. A. M 


. Vermaseren and A. Vogt, Phys. 


625, 245 (2005) 


arXiv:hep-ph/0508055|. 


M. Neubert, 


Eur. Phys. J. C 40, 165 


arXiv:hep-ph/0408179; . 



(2005) 



T. Becher, M. Neubert and B. D. Pecjak, JHEP 0701, 076 
(2007) arXiv:hep-ph/0607228 . 



A. D. 

and R. 



Roberts, W. J. Stirling 

Lett. B 604, 61 (2004) 

arXiv:hep-ph/0410230 . 

L. Magnea and G. Sterman, Phys. Rev. D 42, 4222 (1990). 



Martin, R. G. 
S. Thorne, Phys. 



